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Abstract 

Final results of the search for the Standard Model Higgs boson are presented 
i) . for the data collected by the L3 detector at LEP at centre-of-mass energies up to 

'^ ' about 209 GeV. These data are compared with the expectations of Standard Model 

processes for Higgs boson masses up to 120 GeV. A lower limit on the mass of the 
/\ • Standard Model Higgs boson of 112.0 GeV is set at the 95% confidence level. The 

j^ ■ most significant high mass candidate is a Hz/z/ event. It has a reconstructed Higgs 

mass of 115 GeV and it was recorded at ^/s = 206.4 GeV. 



To be submitted to Phys. Lett. B 



1 Introduction 

One of the most important goals of the L3 experiment at the LEP e"'"e~ colhder was to find the 
Higgs boson. In the Standard Model this particle is associated to the Higgs field jT] , expected to 
provide mass to all the observed elementary particles. The mass of the Higgs boson, mn, is not 
predicted by the theory. Before the advent of LEP, there was no solid experimental information 
about the Higgs mass. The L3 experiment has carried out the search for the Higgs boson at LEP 
in very large data samples collected at the Z resonance [2] and at ever increasing centre-of-mass 
energies and luminosities [3-7] greatly extending the Higgs mass range investigated. 

A fit that includes L3 electroweak precision measurements results in an upper limit on m-a 
of 133 GeV [8^ at the 95% confidence level. Previous L3 direct searches for the Standard Model 
Higgs boson excluded the mass range up to 107 GeV p. Similar results were also reported by 
other LEP experiments ^. Results of the Standard Model Higgs search obtained shortly after 
the end of the LEP data taking in the year 2000 were also published by L3 7J and by the other 
LEP experiments |10j . 

The Standard Model Higgs boson is produced at LEP mainly via the Higgs-strahlung process 
e"'"e~ ^ Z* —>■ HZ. The processes of W^W~ and ZZ fusion contribute, with smaller rate, to 
the Higgs production in the Hi/z/ and He"'"e~ channels, respectively. The largest sources of 
background are four-fermion final states from W and Z pair production, as well as quark pair 
production e"'"e~ ^qq(7). 

In this letter, the final results of the Standard Model Higgs search performed on the data 
collected by L3 at LEP at a centre-of-mass energy, ^/s, up to about 209 GeV are reported. 
These results include the full luminosity collected in the year 2000 and the corrected LEP 
beam energies. In the year 2000 LEP was run at several values of ^/s. The slight beam energy 
adjustments significantly affect the signal expectation at the highest Higgs masses, close to the 
kinematic limit for HZ production, mn = y/s — mz, where mz is the Z boson mass. The effect 
of the Z width in the Higgs mass reconstruction close to the HZ kinematic limit is also taken 
into account. Final calibrations of all subdetectors are applied. The signal and background 
expectations are evaluated on a finer grid of ^/s values, with larger samples of simulated events, 
thus reducing statistical and systematic uncertainties. Therefore, the results reported in this 
letter are affected by total uncertainties smaller than in Reference [7| 

2 Data and Monte Carlo samples 

The data were collected by the L3 detector ^I] at LEP during the year 2000 at several centre- 
of-mass energies. The total collected luminosity amounts to 217.3 pb~^. The data are grouped 
into seven samples corresponding to average centre-of-mass energies between 202.8 GeV and 
208.6 GeV. The integrated luminosities corresponding to these samples are given in Tabled 

The Higgs production cross sections and branching ratios are calculated using the HZHA 
generator J21- Efficiencies are determined using Monte Carlo samples of Higgs events, generated 
with PYTHIA ^3]. Since the Higgs production cross sections and efficiencies depend strongly 
on ^/s, in particular for mn close to the HZ kinematic limit, samples of Higgs events are 
simulated at each centre-of-mass energy shown in Tabled Higgs events are generated with mn 
between 105 and 120 GeV, in steps of 1 GeV. For each mass and each search channel, between 
2000 and 10000 events are generated. 

The Standard Model background estimates rely on the following Monte Carlo programs: 
KK2f ^ for e+e-^qq(7), KORALW [W< for e+e'^W+W", PHOJET [16j for two-photon 



processes (e"'"e~ —>■ e"'"e~qq) and EXCALIBUR 17^ for other four-fermion final states. The 
number of simulated events for the dominant backgrounds is at least 100 times the number of 
events expected for such processes. 

The response of the L3 detector is simulated using the GEANT program P^ , taking into 
account the effects of multiple scattering, energy loss and showering in the detector. Hadronic 
interactions in the detector are modelled using the GHEISHA program ^H]- Time dependent 
detector inefficiencies, as monitored during the data taking period, are also simulated. 

3 Analysis procedures 

The search for the Standard Model Higgs boson is based on the study of four distinct event 
topologies: HZ— i>qqqq, HZ^qqi/z/, HZ — i> qqf+f" (£ = e, /i, r) and HZ-^r+r~qq. 
In the following they are denoted as Hqq, Hz/z/, B.i^i~ and T+r^qq, respectively. With the 
exception of the HZ^r+r^qq decay mode, all the analyses are optimised for the H— i>bb decay. 
This mode represents about 80% of the Higgs branching fraction in the mass range of interest. 

All the search channels are analysed in three stages. First, a high multiplicity hadronic 
event selection is applied to reduce the large background from two-photon processes, while 
preserving most of the Higgs signal. In a second stage, topological and kinematical variables 
together with b-tag variables are either used to construct an event likelihood or fed into a neural 
network, to further discriminate between signal and background events. A b-tag variable is 
calculated for each hadronic jet using a neural network |2j which exploits three-dimensional 
decay lengths, properties of semileptonic b decays and jet-shape variables. The tracking and 
b-tagging performance in the Monte Carlo simulation are tuned using 4 pb~^ of calibration data 
collected at y/s ~ rriz in the year 2000. The b-tagging performance for the high-energy data is 
verified with samples of e"'"e~ —>■ qq(7) events. The efficiency for tagging light flavoured hadrons 
is verifled with W"*"W~ —>■ qqiu events. The agreement of data with the simulation of Standard 
Model processes in the jet b-tag variable based on neural network is shown in Figure Q for 
the e+e~ — >• qq(7) and the W"''W^ -^ qqiu events. The expectation from the Standard Model 
Monte Carlo describes the data within the statistical uncertainty. 

The neural network b-tag variables are combined into an event b-tag variable. First, the 
probability is calculated for each jet to be compatible with the distribution for light quarks 
determined from Monte Carlo. Then, the event b-tag variable is deflned as the negative log- 
likelihood of these probabilities. 

The last part of the analysis is the construction of a flnal discriminant for each topology. 
It is built from a combination of the event likelihood, or the neural network output, with 
the reconstructed Higgs mass. For each Higgs mass hypothesis, the flnal discriminants are 
computed for the data and for the expected background and signal. The distributions of the 
flnal discriminants are then used to calculate the likelihood ratio, Q, as a function of mu. This 
is the ratio of the probability of observing the data in the presence of both the signal and the 
background, "signal+background" hypothesis, to the probability of observing the data in the 
presence of only the background, "background-only" hypothesis. The quantity used to evaluate 
the compatibility of the data with the signal is the log-likelihood ratio deflned by |2nj : 

-2lnQ = 2j2[si - riilnil + Si/bi)] . 

i 

In this expression, i indicates the i-th bin of the flnal discriminant of each channel and at 
each a/s; rij, Sj and bi indicate respectively the number of observed events, the expected Higgs 



signal and the Standard Model background, in the i-th bin. Each event in the sum has a 
weight In (1 + s/b) which depends on the signal-to-background ratio, s/b, in the bin where it is 
found. This weight depends on the Higgs mass hypothesis. For each given mn, the value of the 
log-likelihood ratio in the data is compared to the expected distributions of —2 In Q in a large 
number of simulated experiments under the "background-only" and the "signal+background" 
hypotheses. The results for each search channel are then presented in terms of —2 In Q for the 
data compared to the expected median values for the two hypotheses, as a function of mn. 

3.1 The Hqq analysis 

The Hqq analysis aims to select and study events with four jets, two of which contain b hadrons, 
while the other two must be consistent with the decay of a Z boson. Background from Standard 
Model processes comes mainly from qq final states with hard gluon radiation, W"*"W~ and ZZ 
events, especially those where one of the Z bosons decays into b quarks. 

After a high multiplicity hadronic preselection, the events are forced into four jets with 
the DURHAM algorithm ^21^ and a kinematic fit requiring four-momentum conservation is 
performed. Several discriminating variables, Xj, are combined into a single likelihood which is 
then used to select the final sample. In this combination, each final state is considered as an 
event class j (j=HZ, ZZ, WW, qq). For each class, probability density functions f^{xi) are 
derived from Monte Carlo. The probability for an event to belong to the event class j, based 
on the value of the variable Xi, is defined as p^{xi) = f^{xi)/ ^^ /'^(xj), where k runs over all 
classes. 

The individual probabilities are combined into a likelihood: Lhz = YiiP^'^i^i)/ J2k Yli P^i^i): 
where i runs over all variables considered and k over all event classes. Ten variables are used 
to calculate the likelihood. They are: 

• the number of tracks, 

• the event b-tag, 

• the maximum energy difference between any two jets, 

• the minimum jet energy, 

• the parameter of the DURHAM algorithm for which the event is resolved from three jets 
into four jets, 

• the maximal triple-jet boost, defined as the maximum three-jet boost obtained from the 
four possibilities to construct a one-jet against three-jet configuration in a four-jet event, 

• the minimum opening angle between any two jets, 

• the event sphericity, 

• the mass from a 5C kinematic fit imposing energy and momentum conservation and equal 
dijet masses, M^^ , 

• the absolute value of the cosine of the production polar angle, |cos6|, assuming the 
production of a pair of bosons. 



The distributions of the event b-tag, M^^ , \ cos9| and Lhz for the events selected in the Hqq 
search channel with ^/s > 206 GeV, compared to the expectation for Standard Model processes, 
are shown in Figure El 

Events are selected into the final sample if the value of Lhz exceeds a threshold optimised 
for each centre-of-mass energy and each Higgs mass hypothesis. In addition, the compatibility 
of each event with a Higgs mass hypothesis mu is tested by the variable Xhz = (^« ~ ("^h + 

niz)) /cr|jj2 + {^i — liT^H — fnzl) /"^Ahz' ^^ ^^^^ expression, Ej and Aj are the dijet mass sum and 
dijet mass difference, respectively, for the i-th of the three possible jet pairing combinations, 
while (Tsjj2 ^iid ctahz ^^^ ^^^ corresponding resolutions for Higgs events. The jet pairing with 
the best x^ is chosen. Finally, only events with the x^ probability above 0.01 are selected. As 
an example, for m-n = 110 GeV, 179 events are selected in the data with 172 expected from 
background processes and 12.8 events expected from the Higgs signal; for ttih = 115 GeV, 149 
events are observed with 142 from background and 3.2 from the Higgs signal. 

For these events a final discriminant is constructed. At first, the events are classified into 
three categories depending on the values of the b-tag of the two jets assigned to the Higgs 
boson. The first category contains events where none of these jets has the highest value of the 
b-tag among the four jets of the event. The second category is composed of events where one 
of these jets has the highest b-tag value. The third category contains events where the two jets 
assigned to the Higgs boson have the highest b-tag values. The Xhz probability, the b-tag values 
of the individual jets and the event category are then combined into the final discriminant. 

3.2 The Hvi? analysis 

The Hi/u search is based on the selection of events with two jets containing b hadrons, with large 
missing energy and missing mass consistent with rriz- A neural network is used for the Huu 
analysis, very similar to the one previously reported [HIIZ]. However, tighter cuts on radiative 
photons and on the jet polar angle are applied to reduce the e"'"e~ — *> qq(7) background and 
ensure the best jet energy and b-tag measurements. The signal efficiency is slightly reduced by 
a few percent relative but the search performance is not significantly modified. In addition, the 
neural network is trained with the final, high statistics, signal and background Monte Carlo 
samples, at each ^/s value, to maximise the sensitivity of the analysis. 

In the first step of the analysis, high multiplicity hadronic events are selected and forced into 
two jets using the DURHAM algorithm. The dijet invariant mass must exceed 40 GeV. These 
requirements reduce contributions from two-photon interactions, while retaining a significant 
fraction of hadronic events from e^e" — > qq(7) and W-pair production. These backgrounds 
are then reduced by requiring the visible mass to be less than 140 GeV and the mass recoiling 
against the hadronic system to lie between 50 GeV and 130 GeV. 

Events from e"'"e~ — > qq(7) are further suppressed by requiring the longitudinal missing 
energy to be less than 0.6y^ and the missing momentum vector to be at least 16° away from 
the beam axis. The energy in the forward luminosity calorimeter is required to be below 20 GeV. 
The acoUinearity is required to be smaller than 65°. The distribution of the event b-tag after 
the above cuts is shown in Figure El a. A loose cut requiring the event b-tag to be larger than 
0.5 is then applied, without further loss of signal efficiency. After this set of cuts, there are 123 
events in the data, while 130 are expected from background processes with 4.3 and 1.3 events 
expected for mn = 110 GeV and 115 GeV, respectively. 

A kinematic fit imposing four-momentum conservation and requiring the missing mass to 
be consistent with mz is performed to compute the reconstructed Higgs mass from the two jets. 



The output of a mass independent neural network ^ is then combined with the reconstructed 
Higgs mass to build the final discriminant. The distributions of the reconstructed Higgs mass, 
the missing mass and the neural network output for the events selected in the Huu search 
channel with a/s > 206 GeV, compared to the expectation for Standard Model processes, are 
shown in Figure El General agreement between the data and the expected contributions from 
Standard Model processes is observed in all the distributions. 

3.3 The iii^i~ and T"^T~qq analyses 

The signatures for the He"''e~ and HyU+/i^ processes are a pair of high energy electrons or muons 
with an invariant mass compatible with rriz and two hadronic jets with b quark content. In 
Hr^r^ events the tau pair invariant mass must also be compatible with mz- For these events, 
the mass resolution is worse than in the other B.i~^i~ channels due to the missing neutrinos from 
the tau decays. Events with Higgs decaying into tau leptons, r+r~ qq, have similar signature 
to the Hr+r" events, with the difference that the hadronic jet mass must be compatible with 
mz and that the b-tag content of the event is reduced. 

The analyses are very similar to those described in Reference El The selections require high 
multiplicity events. In the He"'"e~ and H/i^/i" analyses two well identified electrons or muons are 
also required. In the tau analyses, tau leptons are identified either by their decay into electrons 
or muons, or as an isolated low-multiplicity jet with one or three tracks and unit charge. The 
identified leptons must have a large opening angle and must be well isolated from the hadronic 
jets. For all H£"^£~ selections, the invariant mass of the leptons after a kinematic fit imposing 
four-momentum conservation must be consistent with mz within a mass range depending on 
the dilepton mass resolution. In the r"'"r~qq selection the mass of the two hadronic jets after 
kinematic fit must be consistent with m^z- 

After the B.i~^i~ selection, 18 events are observed with 16.7 expected from background 
processes and 1.7 or 0.32 signal events expected for mn =110 GeV or 115 GeV, respectively. 
After the r+r'qq selection, 8 events are observed with 7.8 expected from background and 0.66 
or 0.15 signal events expected for mu =110 GeV or 115 GeV, respectively. 

The distributions of the dilepton mass and the reconstructed Higgs mass in the He"'"e~ and 
H/x"*"/!" channels are shown in Figure |31i and|3jD. The distributions of the reconstructed Higgs 
mass in the Hr+r^ and r+r^qq channels are shown in Figure ID:; andlUi, respectively. 

In the B.i~^i~ selections, the dijet mass after the fit is combined with the b-tag values of the 
two jets, to form the final discriminant. For the r"'"r~qq selection, the mass of the tau pair, 
calculated by constraining the invariant mass of the two other jets to mz, is used as the final 
discriminant. 

4 Results 

Figure El shows the observed —2 In Q compared to the expectation for the "background-only" 
and the "signal+background" hypotheses, as a function of mn, for each of the search channels. 
An observed value of —2 In Q larger than the median expected value for the background indicates 
a deficit of events with respect to the expected background while an observed —2 In Q value 
below the median expected background value indicates an excess. Good agreement between 
the observation and the expected background is observed in all channels within one standard 
deviation from the background expectation. A slight excess of events above one standard 
deviation from the background is observed in the Hi/z/ channel for mn above 100 GeV. The 



observed and expected log-likelihood ratio —2 In Q for all channels combined as a function of 
rriH is shown in Figure IHl 

These results are used to evaluate confidence levels for the "background-only" and the 
"signal+background" hypotheses. The confidence level for the "background-only" hypothesis, 
1 — CLft ^20j, is the probability of observing in a large sample of simulated "background- 
only" experiments a more signal-like value of the log-likelihood ratio than is actually observed. 
The distribution of 1 — CL;, in a large sample of "background-only" experiments is uniform 
between and 1, thus its median expected value is 0.5. An observed value of 1 — CL;, 
lower than 0.5 indicates an excess of events in data compared to the expected background. 
Similarly, the "signal+background" confidence level CL^+b is defined as the probability in a 
sample of "signal+background" experiments of observing a less signal-like value of the log- 
likelihood ratio than is actually observed. To exclude a signal, an additional quantity is defined, 
CLs=CLs+b/CLb [211]. The signal hypothesis is excluded at 95% confidence level when CL^ has 
a value smaller or equal to 5%. 

The statistical and systematic uncertainties on the signal and background expectations are 
included in the calculations of the combined confidence levels. Statistical uncertainties on the 
background and signal predictions, arising from the finite number of generated Monte Carlo 
events, are evaluated to be up to 8% for the background and 4% for the signal. The systematic 
uncertainties are derived using a similar procedure to the one adopted in previous Standard 
Model Higgs searches 0. In addition, a systematic uncertainty on the qq background, which 
affects mostly the search region close to the HZ kinematic limit in the Yiuv and Hqq channels, 
is included depending on mn- Thus the systematic uncertainty on the number of background 
events is estimated to be from 6% up to 15% for rran close and beyond the HZ kinematic limit. 
The systematic uncertainty on the number of signal events is estimated to be between 3% and 
6%, for mn close to and beyond the HZ kinematic limit, to take into account the spread of ^/s 
values in the different data samples. 

The statistical uncertainty is uncorrelated from bin to bin in the final discriminant distri- 
butions and has little effect on the confidence level. Bins of the final discriminant distributions 
with a s/h ratio below 0.05 are not considered in the calculation of the confidence levels, as they 
degrade the search sensitivity once systematic uncertainties are included in the calculation. The 
number of selected and expected events for all the analyses after such a s/h cut are summarised 
in Table Elfor the data, the background and the Higgs signals for mn = HO GeV and 115 GeV. 
The number of signal events includes cross-efficiencies from other channels, fusion processes 
and charm and gluonic Higgs decays. 

The confidence level for the "background-only" hypothesis 1-CLb and the confidence level 
for the signal hypothesis CL^ as a function of mn are shown in Figure [7| They are computed 
following the procedure of Reference 1201 The results of the L3 Standard Model Higgs searches 
at lower centre-of-mass energies [3101 are included in the calculation of these confidence levels. 
Values of mn below 107 GeV are excluded in the Standard Model with a confidence level greater 
than 99.5%. 

The observed lower limit on tt^h is 112.0 GeV at the 95% confidence level, for an expected 
lower limit of 112.4 GeV. This new value improves upon and supersedes our previously pub- 
lished limit. For mn =112.0 GeV, where CL^ is 5%, the background probability l-CLf, is 40%. 
For mn =115 GeV, the background probability is 32%. The previously published background 
probability estimates [7| are consistent with the final results presented here, given the size of 
the uncertainties affecting the signal and background estimate in the vicinity of the kinematic 
limit. 



The most significant candidate for mn =115 GeV is a Huu event. It lias a reconstructed 
Higgs mass of 115 GeV and it was recorded at ^/s =206.4 GeV. The kinematic properties of 
this event were described in detail in Reference El 
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i (GeV) 
Luminosity (pb" 



202.8 203.8 205.1 206.3 206.6 208.0 208.6 
2.7 7.6 68.1 66.9 63.7 8.2 0.1 



Table 1: The average centre-of-mass energies and the corresponding integrated luminosities of 
the data samples collected in the year 2000. 



200< v^ <209 GeV 


Mass hypothesis 


mu 


= 110 GeV 


mH = 115 GeV 


Selection 


Nd 


Nb Ns 


Nd Nb Ns 


Hqq 


49 


51.5 11.7 


12 9.4 1.8 


Hz/z/ 


13 


10.7 3.3 


5 3.3 0.66 


He+e- 





0.66 0.58 


0.38 0.14 


Rfi+fi- 





0.38 0.45 


0.26 0.11 


Hr+r- 





0.53 0.19 


1 0.14 0.03 


r+r- qq 


3 


2.3 0.51 


0.84 0.15 


Total 


65 


66.1 16.7 


18 14.3 2.9 



Table 2: The number of observed candidates (Nd), expected background (Nb) and expected 
signal (Ns) events for the data collected in the year 2000, after a cut on the final discriminant 
corresponding to a signal-to-background ratio greater than 0.05. This cut is used to calculate 
the confidence levels. 
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Figure 1: Distribution of the neural network jet b-tag variable in a sample of a) e^e~ — >qq(7) 
and b) W^W~ -^ qqiu events selected from the high-energy data collected in the year 2000. 
Two entries per event contribute to the distributions. The data are compared to the simulation 
of Standard Model processes. The bin-by-bin ratio of the data to the simulated events is 
displayed in c) and d). 
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Figure 2: Distribution of a) the event b-tag, b) the cosine of the boson production angle, c) 
the mass from the 5C equal-mass fit and d) the likelihood for the events selected in the Hqq 
search channel. The points correspond to the data collected at y/s > 206 GeV. The open and 
hatched histograms are the expected backgrounds from Standard Model processes. The dashed 
line is the distribution expected for a 115 GeV Higgs signal, multiplied by a factor of 200. 
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Figure 3: Distribution of a) the event b-tag, b) the reconstructed Higgs mass, c) the missing 
mass and d) the neural network output, for the events selected in the Hz/z/ search channel. The 
points represent the data collected at ^/s > 206 GeV. The open and hatched histograms are 
the expected backgrounds. The dashed line is the expected Higgs signal with mu = 115 GeV, 
multiplied by a factor of 30. 
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Figure 4: Distributions of a) the dilepton mass and the reconstructed Higgs mass in the b) 
He"'"e~ and H/i"'"/i~, c) Hr"'"r~ and d) r"'"r~qq channels. The points are the data and the open 
and hatched histograms the expected backgrounds. The dashed line is the expected Higgs 
signal with mu = 115 GeV, multiplied by a factor of 30, in each channel. 
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Figure 5: The log-likelihood ratio, —2 In Q, as a function of the Higgs mass hypothesis, mn, for 
the search channels a) Hqq, b) Huu, c) B.i~^i~ and d) r"^r~ qq. The solid line shows the observed 
—2lnQ. The dashed line shows the expected median value of — 21nQ for the "background- 
only" hypothesis. The dark and light shaded bands show the 68% and 95% probability intervals 
centred on the background expected median value. The dotted line is the median expected value 
for the "signal+background" hj^othesis. 
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Figure 6: The log-likelihood ratio, —2\nQ, as a function of the Higgs mass hypothesis, mn, 
for all the search channels combined. The solid line shows the observed — 21nQ. The dashed 
line shows the expected median value of —2lnQ for the "background-only" hypothesis. The 
dark and light shaded bands show the 68% and 95% probability intervals centred on the back- 
ground expected median value. The dotted hne is the median expected value for the "sig- 
nal+background" hypothesis. 
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CLfc and b) the signal confidence level, CL^, 
for all the search channels combined. The 



Figure 7: a) The background confidence level 1 - 
as a function of the Higgs mass hypothesis, mn 
data collected at 189 < ^/s < 202 GeV |,5ill| are also included in the combination. The solid 
line shows the observed value. The dashed line shows the median expected value in a large 
number of simulated "background-only" experiments. The dark and light shaded bands show 
the expected 68% and 95% probability intervals centred on the background expected median 
value. The observed lower limit on the Higgs mass is set at 112.0 GeV, with an expected median 
value of 112.4 GeV, at the 95% confidence level. 
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